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We report precision measurements of the effective mass m* in high-quality bilayer graphene using 
the temperature dependence of the Shubnikov-de Haas oscillations. In the density range of 0.7 
x 10 12 /cm 2 < n < 4.1 x 10 12 /cm 2 , both the hole mass and the electron mass mt increase 
with increasing density, demonstrating the hyperbolic nature of the bands. The hole mass 
is approximately 20-30% larger than the electron mass m*. Tight-binding calculations provide a 
good description of the electron-hole asymmetry and yield an accurate measure of the inter-layer 
hopping parameter V4 = 0.063. Both m£ and m* are suppressed compared to single-particle values, 
suggesting electron-electron interaction induced renormalization of the band structure of bilayer 
graphene. 

PACS numbers: 73.43. Qt, 73.20.At, 71.70.Gm. 



I. INTRODUCTION 

Bilayer graphene may be a technologically important 
material in electronics and photonics due to its tunable 
band gap. The fundamental property that underpins 
such applications-its band structure-has been the sub- 
ject of many recent theoretical^ and experimental stud- 
ies using angle-resolved photoemission spectroscopy^, in- 
frared and Raman measureme nts^^ , cyclotron mass^ and 
compressibility measurements^^. On a single parti- 
cle level, the band structure of bilayer is thought to be 
well described by a tight-binding HamiltoniarP^l w ith 
a few leading order Slonczewski-Weiss-McClure param- 
eters, i.e., 7o, 71, 73 and 74. Experimental knowl- 
edge of these hopping parameters in bila yer varies, with 
71 = 0.40 eV fairly accurately knowrF^ anc j t ne res t to 
a much less degree. For example, experimental values 
of 74, which cont rols the band asymmetry, range from 
O.H-0.19 eVEEED 

Meanwhile, electron-electron interactions in single- 
layer and bilayer graphene are predicted to be strong 
and peculiar. Inte resting collective states emerge in a 
magnetic fiel d 12 * 13 -! The many-body corrections to Fermi 
liquid parameters such as the compressibility k and the 
effective mass to* are expected to be sunstantial already 
at currently accessible densitie^HE These renormaliza- 
tion effects are related to, but also quantitatively differ- 
ent from those observe d in co nventional two-dimensional 
electron gases (2DEGs) 19 20 , due to the chirality of sin- 
gle and bilayer graphene^. For example, instead of an 
enhancement!^, the effective mass of bilayer graphene is 
predicted to be increasingly suppressed at lower carrier 
densities^. No experimental evidence of such renormal- 
ization effect has been reported so far. 

In this work, we report the measurements of the effec- 
tive mass m* in bilayer graphene samples for a wide range 
of carrier densities using high-quality Shubnikov-de Haas 
(SdH) oscillations. The inter-layer hopping parameter 



74 is determined to be 74 = 0.063(1)7 with the highest 
accuracy reported so far. The magnitude and density de- 
pendence of to* deviate from tight-bind calculations, pro- 
viding evidence for electron-electron interaction induced 
band renormalization. 



II. SAMPLE PREPARATION 

Bilayer graphene flakes are exfoliated onto 290 nm 
Si02/Si wafers from highly ordered pyrolytic graphite 
and identified by optical microscope and Raman spectra. 
They are further confirmed by their quantum Hall se- 
quence. Conventional electron-beam lithography is used 
to pattern the flakes into Hall bars. 



III. MEASUREMENTS 

Transport measurements are carried out in a He 4 sys- 
tem using standard low-frequency lock-in technique. The 
field effect mobility /^fe = (l/e)(d<r/dn) of our pristine 
bilayer graphene ranges from 3 000-12 000 cm 2 /Vs. Data 
from two samples (A and B) are presented in this paper. 

In Fig. [TJ we plot the sheet conductance a vs. the 
back-gate voltage Vb g of sample A at selected tempera- 
tures between 15 K-250 K. At 15 K, the mobility /xfe 
of sample A is approximately 4 800 cm 2 /Vs for holes 
and 3 100 cm 2 /Vs for electrons. Sample B has a higher 
mobility of 6 300 cm 2 /Vs for holes and 6 800 cm 2 /Vs 
for electrons. The conductance of bilayer graphene sam- 
ples shows a variety of temperature dependence, depend- 
ing on the carrier density and mobility. Near the charge 
neutrality point, all our samples show an insulating-like 
T-dependence (dcr/dT > 0), as shown in Fig. [I] This 
behavior is due to the thermal excitation of carriers out 
of electron-hole puddles, as demonstrated in Ref. I2T1 As 
the carrier density increases, dcr/dT eventually becomes 
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FIG. 1. (color online) (a) and (b) The sheet conductance 
cr(Vbg) of sample A. From top to bottom: T= 15, 100, 150, 
250 K in (a) and in reverse order in (b). The charge neutrality 
point is at Vb g = 7 V. 

negative (metallic) in the highest-quality samples. This 
trend is shown by the hole branch in Fig. [lja) and (b), 
where the crossover density is approximately = 2.1 
x 10 12 /cm 2 . For samples with lower mobilities, the 
insulating-like T-dependence persists to high densities, 
an example of which is given by the electron branch in 

Fig.njb). 

This complex behavior is in contrast to that of single- 
layer graphene, where a metallic-like temperature depen- 
dence dominates over a wide range of densities due to 
phonon scatteringpHU! qualitative features of our 

data are consistent with the model proposed in Ref. [SSI 
where cr(T) combines metallic and insulating trends aris- 
ing from the conduction of the majority and minority 
carriers respectively. The true metallic T-dependence of 
a bilayer graphene two-dimensional electron gas emerges 
only in high quality samples and/or at high carrier den- 
sities. In Fig. [I] the different T-dependence of the two 
carriers in the same sample points to an intrinsic electron- 
hole (e-h) asymmetry of bilayer graphene, which we fur- 
ther examine below. 

To probe the band structure of bilayer graphene, we 
measure the effective mass m* as a function of the carrier 
density using SdH oscillations. This technique is well es- 
tablished in 2DEGs but require high-quality oscillations 
to reliably extract m* . Figure [2ja) shows the SdH oscil- 
lations p xx (i?) of sample A at a high electron density n e 
= 3.26 x 10 12 /cm 2 and varying temperatures. The oscil- 
lations have an early onset, appear sinusoidal and free of 
beating. Its amplitude 5p xx is given by-23 

6/Q xx 7r 27r 2 k B T/hui c 

= 4 7 t h exp( ); 7 t h = . ,, n 2 . , (1) 

po u> c T q smn(27r z KBJ /nu! c ) 

where u c = eB/m* is the cyclotron frequency, r q is the 
quantum scattering time and "f t h the thermal factor. 

As shown in Fig. [2^a) , <5p xx increases with increasing 
B and decreasing T . Its T-dependence provides a direct 
measure of m* whereas the independence is controlled 
by both m* and r q . At each carrier density, the low-field 
<5p xx (T, B) data, i.e., before the onset of the quantum 
Hall effect, is fit to Eq. [Ijwith two fitting parameters m* 
and r q . The simultaneous fitting of m* and r q allows 
us to accurately determine <5p xx , especially at low car- 
rier densities, where the oscillations are few and a linear 
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FIG. 2. (color online) (a) SdH oscillations p xx (-B) at T = 1.5- 
50 K for n c = 3.26xl0 12 /cm 2 . (b) p xx (B) at T = 10 K and 40 
K for rih = 0.89 x 10 12 /cm 2 . Dashed lines are fittings with r q 
= 42 fs and ra£ = 0.036m c . A smooth background has been 
subtracted, (c) <5p xx /T vs. T in a semi-log plot for n c — 3.26 
x 10 12 /cm 2 at B = 7.53 T (down triangle in (a)) and for rtj, = 

0. 89xl0 12 /cm 2 at B = 4.70 T (up triangle in (b)). The sym- 
bols correlate. Dashed lines are fittings with m* = 0.041m o 
(down triangles) and m h = 0.036m c (up triangles). (d) r q vs. 
n for electrons (red triangles) and holes (black squares). All 
data in (a)-(d) are from sample A. 

interpolation between peaks, as commonly done in the 
literature^!, cannot give the correct amplitude of <5p xx . 
Figure [2jb) shows p xx (S) data at T=10 K and 40 K for 
a low hole density = 0.89 x 10 12 /cm 2 . Fittings to Eq.[l] 
are shown as dashed lines. Only the right values of m* 
and T q can fit both the independence and T-dependcncc 
of <5p xx simultaneously. 

In Fig. [2]^c), we plot two examples of the measured 
Spxx/T vs. T in a semi-log plot for the two positions 
marked in Figs. [2]ja) and (b) with down and up triangles 
respectively. Dashed lines are fittings generated with m* 
= 0.041to c and m h = 0.036m c respectively, where m is 
the electron rest mass. They both fit very well. Overall, 
Eq. [l] provides an excellent description of the Sp xx (T, B) 
data in the entire density range studied, with the un- 
certainty of rn* increasing from 0.0001to c to 0.0015rn c 
from high to low densities. The global fitting proce- 
dure also ensures that the extracted m* is filling factor- 
independent and therefore represents the B = limit, 

1. e., the band structure mass. This m* is not directly 
comparable to m* determined from cyclotron resonance 
measurements^, as Coulomb interaction may manifest 
differently in these two cases^SI. A good illustration of 
this situation is the parabolic band material GaAs, where 
m* determined from SdH oscillations embodies electron- 
electron (e-e) interaction 2 ^ while its effect is forbidden 
in cyclotron resonance measurements by the Kohn theo- 
reirPSl. 

Using this method, we have determined m* and r q for 
samples A and B in the density range of 0.7 < n < 4.1 
x 10 12 /cm 2 for electrons and holes. Both samples show 
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oscillations of equally high quality and comparable r q 
. Figure [21(d) plots r q (n) of sample A for both carri- 
ers. Overall, r q increases with increasing density, ranging 
from 41 to 60 fs. These values correspond to a disorder 
broadening T = h/2T q = 5.5-8.0 meV, which are si mila r 
to high-quality single and bilayer graphene sample^^fil. 

The results of m* of samples A and B as a function of 
the carrier density n are plotted in Fig.ga). The error 
bars represent uncertainties obtained from fittings simi- 
lar to those shown in Fig. |2jc). The two samples agree 
very well with each other. In the density range studied, 
both m* and increase with increasing n, indicating 
the non-parabolic nature of the bands. This observation 
agrees with the compressibility measurements of Refs.HOl 
and [11] and is also consistent with the observation of a 
constant m* at yet lower densities^. The ratio of ra^/m* 
is about 1.2-1.3, demonstrating a pronounced electron- 
hole asymmetry. 



IV. ANALYSIS AND DISCUSSION 

The above measurements of m* provide an accu- 
rate means of determining the band structure of bi- 
layer graphene and investigating the effect of e-e interac- 
tion. In the following analysis, we employ a tight-binding 
Hamiltonian following the notations of Refs. |3]and[7J 

/V(n)/2 + A 4> 7i ~v4(f>* \ 

_ | <f>* V(n)/2 -v 4 <f>* v 3 

7i -v±cj> -V(n)/2 + A <p* 

V -v A <f> v 3 <t>* </> -V(n)/2J 

(2) 

Equation[5]is written in the basis of the four sublattices 
(^Ai: ^bi, *A2, *B2)j where Al, A2 are the two stacked 
sublattices in layer 1 and 2 respectively. The nearest 
neighbor in-plane (Al-Bl) hopping integral 70 is included 
in (j) — 7o(3/2/c y a — i3/2fc x a) = hvp(k y — ik x ), where a 
= 1.42A is the carbon-carbon distance and the momen- 
tum vector (k x ,k y ) originates from the K (K 7 ) point of 
the Brillouin zone. The Fermi velocity vp = (3/2)joa/h. 
7i, v 3 = 73/70 and v-i — 74/70 represent the hopping 
integrals between two inter-layer sublattices A1-A2, Bl- 
B2, and A1-B2 respectively. 71 gives rise to the band 
splitting, 73 leads to trigonal warping of the Fermi sur- 
face and 74 controls electron-hole asymmetry. A is the 
on-site energy difference of Al and Bl, due to their stack- 
ing difference. V(n) is the potential difference between 
the two layers and varies with the carrier density^. The 
eigenvalues of Eq.[2]produce the four low-energy bands of 
bilayer graphene. Out of the four bands, the two higher 
energy electron and hole bands are neglected here since 
they are far above the Fermi level of our density range, 
E-p <~ 30-120 meV. The effective mass m* of the lower 
bands is given by: 

* ^ dA(E) 



TABLE I. The effect of tight-binding parameters on m* and 
their values. 
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where A(E) is the fc-space area enclosed by the contour 
of constant energy E. For 73 = 0, the contour is circular 
and Eq.[|is simplified to m* = h 2 k / (dE (k) / dk) . 

We diagonalize Eq. [2] and numerically compute m* us- 
ing Eq. [3j The effect of each parameter in Eq. [2] on m* is 
summarized in table [TJ where +(-) means an increase of 
the parameter will increase (decrease) the value of m* . 

In our calculations, the inter-layer B1-B2 hopping en- 
ergy 73 is set to zero due to its negligible effect in the 
density range considered here (see Appendix A for de- 
tails). The gate voltage-induced V(n) is calculated fol- 
lowing Eqs.(7)-(13) in Ref.[?l using self-consistent screen- 
ing and including the small initial doping of our samples. 
Both V(n) and the initial doping produce minute correc- 
tions to m* in the density range studied (sec Appendix B 
for details). Consequently, the overall magnitude of m£ 
and 77i* and their density dependence are predominantly 
controlled by 70 and 71. In the literature, 71 is found 
to be 0.38-0.40 eV by infrared measurement^^. Most 
of our fittings use 71 =0.40 eV. Alternative scenarios are 
also considered in the discussion of electron-electron in- 
teraction effect and further explored in Appendix C. The 
difference between and m* is controlled by and A. 
We fix A=0.018 eV in our calculations. A 10% variation 
of A among literature values^' leads to a change smaller 
than 2% in i>4, which is comparable to its estimated un- 
certainty. 

Fitting to 771^ and m* simultaneously allows us to de- 
termine the remaining adjustable parameters, 70 and V4. 
Fittings to both samples are given in Fig. U a). The value 
of 70 varies slightly from 3.419 eV in A to 3.447 eV in 
B, yielding an average 7o=(3.43 ± 0.01) eV. This cor- 
responds to a Vp = 1.11 x 10 6 m/s, in agreement with 
previous experiments. Both samples yield i>4 = 0.063 ± 
0.001. The value of V4 is also independent of the choice 
of 71, as the fittings in Fig. |3][a) show. This result is 
consistent with the range of V4 ^0.04-0.06 obtained pre- 
viouslySHgjuy k u f. nag a muc h higher precision. This ac- 
curate knowledge of electron-hole asymmetry will be im- 
portant to potential electronic and optical applications 
of bilayer graphene. 

The above fitting does not include the in-plane next- 
nearest neighbor hopping integral 7 n ^, which also con- 
tributes to the electron-hole asymmetry of m* , acting in 
the opposite direction of vj*^. The value of 7„ is not 
well established. Including an additional diagonal term 
— 7n I 4> | 2 /7i in Eq. [2p our calculations show that the 
effect of 7 n on ?j 4 can be represented by an empirical re- 
lation 7J4 = 0.063 + 0.0377 n , which can provide further 
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correction to should the value of 7 n become known. 

The fittings in Fig. [3^ a) reveal an important trend of 
our data, i.e., the measured m* increasingly drops below 
the calculated m* as n decreases. This trend is consis- 
tently seen for both electrons and holes and in both sam- 
ples. Extensive tests show that this discrepancy between 
data and tight-binding calculations cannot be reconciled 
by varying any other parameters except for 71 . A perfect 
fit to both high and low densities is only possible if 71 
is allowed to decrease from 0.40 eV to 0.30 eV, as shown 
by the short-dashed lines in Fig. |3]ja). This scenario, 
although appealing, is at odds with previous experimen- 
tal determination of 7i=0.38-0.40 eV from infrared spec- 
troscopy^^. Alternatively, we attribute the discrepancy 
between measurements and tight-binding calculations to 
interaction induced band renormalization effect. Indeed, 
a recent calculation in bilayer graphene predicts a mono- 
tonic suppression of m* as a function of decreasing den- 
sity^ and the effect is shown to be already substantial 
in the density range studied here. First-principles calcu- 
lations also show that a more complete inclusion of e-e 
interaction in the form of GW corrections increases 70 
from the mean-field-like value o f 2.7 eV to an interaction- 
modified value of 3.4 e V 4 * 31 * 32 l In our experiment, the 
suppression of m*, its density dependence, and the fit- 
ting result of 7o=3.43 eV all strongly point to the renor- 
malization effect of e-e interaction on m* . 

The magnitude of this effect is illustrated in Fig. [3jb) 
and (c), using sample A as an example. Here we calculate 
and plot three sets of values using 70 = 2.72, 3.09, 
and 3.42 eV (corresponding to v F = 0.88, 1.0, 1.11 x 10 6 
m/s respectively). These three values are the first prin- 
ciple mean-field-like, intermediate, and our fitting result 
of 70 respectively. The other parameters are fixed at 
values listed in Table [i] Fig. |3jc) plots the ratio of mea- 
sured m* and the calculated TOq, m* jm§ vs n for each 70. 
The trend of decreasing to*/toq with decreasing density 
is seen in each case, with the magnitude of the suppres- 
sion depending on the input value of 70. Electrons and 
holes follow the same trend. When the first principle 
mean-ficld-likc value of 70 = 2.72 eV is used (bottom 
trace in Fig. [3jc)), the suppression of m* is quite large, 
varying from 0.6 to 0.7 in the density range 0.7 < n < 4 
x 10 12 /cm 2 . These observations provide the first experi- 
mental indication of e-e interaction induced band renor- 
malization effect in bilayer graphene. The quantitative 
input provided by our data should constrain and guide 
future calculations on this important subject, as the cor- 
rect theory must capture both the magnitude and the 
density dependence of m* . 



V. CONCLUSION 

To conclude, we report the measurement of the ef- 
fective mass m* in bilayer graphene over a wide range 
of electron and hole densities. Our results demonstrate 
a pronounced electron-hole asymmetry, from which we 
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FIG. 3. (color online) (a) Measured mj^ and ml vs. n for sam- 
ples A and B. The symbols are indicated in (a) and used in 
(a)-(c). Solid blue lines are fittings to sample A with 70 = 
3.419 eV, 71 = 0.40 eV and u 4 =0.063. The magenta dash- 
dotted line is a fitting to sample B with 70 = 3.447 eV, 71 = 
0.40 eV and U4=0.063. The yellow short dashed lines corre- 
spond to 70 = 3.167 eV, 71 = 0.30 eV and u 4 =0.063. (b) A 
comparison of measured m* and calculated mj for sample A. 
From top to bottom: 70 = 2.72 (olive), 3.09 (wine), and 3.42 
eV(blue). 71 = 0.40 eV, u 4 =0.063 and A = 0.018 eV for all 
traces, (c) The ratio m* /mo vs. n for sample A. From top to 
bottom: 70 = 3.42, 3.09, and 2.72 eV. Dashed lines are guide 
to the eye. 

accurately determine the inter-layer hopping parameter 
V4 in the tight-binding description of the band struc- 
ture. The measured m* is suppressed compared to 
single-particle predictions, indicating the possibility of 
interaction-induced band renormalization at play. Our 
results provide critical experimental input to understand 
the effect of electron-electron interaction in this unique 
two-dimensional electron system. 
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Appendix A: THE EFFECT OF 73 ON m* 

The inter-layer B1-B2 hopping energy 73 leads to trigo- 
nal warping of the Fermi surface^. An example is plotted 
in Fig. ||| for n ~ 0.7 x 10 12 /cm 2 and n ~ 4 x 10 12 /cm 2 . 
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FIG. 4. (color online) The warped Fermi surfaces in momen- 
tum space for E F = 100 meV (n~4x 10 12 /cm 2 ) (red line) 
and E F = 25 meV (n ~ 0.7 x 10 12 /cm 2 ) (blue line). The 
parameters used are 70 = 3.1 eV, 71 = 0.4 eV, 73 = 0.31 eV, 
A = and V(n) = 0. 



Appendix B: THE EFFECT OF V{n) AND THE 
INITIAL DOPING ON m* 

The gate voltage-induced V(n) is calculated using 
Eqs.(7)-(13) in Ref. [71 which agrees well with optical 
measurements^!. The calculation includes a small ini- 
tial hole doping of Vb g = +7 V for sample A and Vb g = 
+17 V for sample B, likely due to water adsorbed on the 
top layer. It also includes the quantum level broadening 
r calculated from the quantum scattering time r q . V(n) 
increases with increasing n and varies from 2-70 meV. As 
V(n) increases, both m£ and to* increase slightly com- 
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FIG. 5. (color online) The measured (symbols) and calculated 
(solid lines) m* of sample A. Blue lines correspond to an initial 
doping of Vbg = +7 V and orange lines correspond to zero 
doping. The other parameters used are 70 = 3.419 eV, 71 = 
0.4 eV, 73 = 0.0 eV, v 4 = 0.063 and A = 0.018 eV. 



to* is calculated using Eq. [3] of the main text. In the 
density range of our experim ent, e ven the largest 73 = 
0.38 eV found in the literatur e 1 5 * 6 * 8 ! results only a minute 
increase of to* , comparable to the smallest error bar. We 
therefore set 73 = in the calculations. 



FIG. 6. (color online) Fitting results of 70 as a function of 71. 
The linear fit corresponds to 70 = 2.411 + 2.5271 . 

pared to zero-gap cases. The effect is the largest at V 
= 70 meV (n=4.1 x 10 12 /cm 2 ), where both to^ and to* 
increase by ~ 0.002 m e , which is ~ 5 % of measured to*. 

Fig. [5] illustrates the effect of the initial chemical dop- 
ing on the calculated to* in sample A. The blue lines 
are calculated with the measured doping of Vb g = +7 V. 
The orange lines are calculated with zero doping. The 
initial doping drastically enhances to* near the charge 
neutrality point but produces negligible effect in the den- 
sity range studied here. 



Appendix C: THE CORRELATION BETWEEN 70 
AND 71 

In Fig. |3ja) , the fitting results of 70 depend on the in- 
put parameter ji . This relationship can be described by 
a linear fit as shown in Fig. [6] As discussed in the text, 
the choice of 71=0. 40 eV leads to 7o=3.419 eV and dis- 
crepancy between data and calculations at low densities. 
A decrease in 7! also decreases 70 and leads to better fit 
at low densities. Both high and low density data can be 
fit by 7i=0.30 eV and 7o=3.167 eV. This choice of 71 
is however incompatible with the experimental range of 
0.38-0.40 eV obtained from infrared absorption measure- 
ments^. 
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